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Abstract 

Novel analytic solutions are derived for integrals that involve the generalized Marcum Q—function, 
exponential functions and arbitrary powers. Simple closed-form expressions are also derived for spe¬ 
cific cases of the generic integrals. The offered expressions are both convenient and versatile, which 
is particularly useful in applications relating to natural sciences and engineering, including wireless 
communications and signal processing. To this end, they are employed in the derivation of the average 
probability of detection in energy detection of unknown signals over multipath fading channels as well 
as of the channel capacity with fixed rate and channel inversion in the case of correlated multipath 
fading and switched diversity. 
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I. Introduction 


The generalized Mareum Q—funetion, Qm{ci,b), has been extensively involved in numer¬ 
ous areas of wireless communieations including digital communications over fading channels, 
information-theoretic analysis of multi-antenna systems, cognitive radio communications, radar 
systems, m-m and references therein. Furthermore, its use has enabled the derivation of several 
tractable analytic expressions for important performance measures in communication theory [|4l|. 

The derivation of tractable analytic expressions in natural sciences and engineering is typically 
a tedious, if not impossible, task because cumbersome integrals are often encountered ll9]l- 
[l22]l . This is also the case when the Marcum Q—function is involved in integrands along with 
exponential functions and arbitrary power terms. Two such integrals are the following: 



( 1 ) 


and 



( 2 ) 


These integrals have been widely employed in the analysis of multi-channel receivers with non¬ 
coherent and differentially coherent detection as well as in the detection of unknown signals 
in cognitive radio and radar systems [Q, Il2^ - ll3^ and the references therein. Based on this, a 
recursive formula for (l2l), that is restricted to only integer values of k and m, was firstly reported 
in dH. Likewise, exact infinite series for ([T]) and © were proposed in ifT^ while a closed-form 
solution to @ for integer values of k was recently reported in ll2T1l . 

Nevertheless, the existing expressions for O and © are subject to validity restrictions, which 
limit the generality of the involved parameters and often render them inconvenient to use in 
applications of interest. Motivated by this, the present work is devoted to the derivation of 
novel closed-form expressions for © and ©, which are more generic and have a relatively 
tractable algebraic representation. These characteristics render them useful in several analyses 
in natural sciences and engineering, including the broad areas of wireless communications and 
signal processing. To this end, they are subsequently employed in the derivation of closed-form 
expressions for the following indicative applications: i) the average probability of detection in 
energy detection over Nakagami—m multipath fading channels which, unlike previous analyses. 
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is valid for arbitrary values of m; ii) the channel capacity with channel inversion and fixed 
rate in arbitrarily correlated Nakagami—m fading conditions using switch-and-stay combining. 
The derived expressions are validated extensively through comparisons with respective computer 
simulations results. 

II. Analytic Solutions to Integrals Involving Power, Exponential and 

Marcum Q-functions 

A. Closed-form Solutions to Q{k, m, a, b,p) 

Theorem 1. For a, 6 G M, m G N and k,p ^ M’*', the following closed-form representation holds 



(3) 



where r(-), and iFif) denote the Euler Gamma function, the Humbert hypergeometric 

function of the first kind and the Kummer hypergeometric function, respectively / |3^ . 

Proof: The series in [fT^ eq. (10)] can be re-written as followi] 


g{k,m,a,b,p) = ^ 


m—1 


T{k -\- n) iFi ( A; + n; n + 1 






(4) 


A{k,a,b,p) 


Using the infinite series in [|4^ eq. (9.14.1)], it follows that 


It —u t=u 

where (x)„ denotes the Pochhammer symbol. Given that 



(5) 


'Equation (10) in QU contains a typo since the (a^ + 2p) term in the denominator should read as [b^ + 2p). This has been 
corrected in 
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and 


one obtains 


{k + n)i 


{k^n+l 

{k)n 




(l)n+i 

(l)n 


A{k,a,b,p) = 

n=0 1=0 


{k)n+lA)n \b'2+‘^P J \2f)2+4p 
{l)n+i n\ l\ 


where 

2^r(fc)exp(aV2) 

(62 + 2pf 

Notably, eq. ([8]) ean be expressed in terms of the Humbert funetion, $i, namely 


( 6 ) 

(7) 


( 8 ) 

(9) 


A(k,a,b,p) = <>■ 1;^.^^) (10) 

Inserting (fTOl) in dH) yields Q, whieh eompletes the proof. ■ 

It is noted that Humbert funetions and their properties have been studied extensively over the 
past deeades ll^ . Il40ll . 


Theorem 2. For a, 6 G M and m,p & and k eN, the following closed-form expression holds 


Q{k,m,a,b,p) 


T{k) T{k)b^^e °2 ^ (m);(2p)^ / \ 

pk p^(62 + 2p)”^/!(62 + 2p)* ^ \ ’ ’ 262 + 4p/ 


Proof: By integrating dB by parts, it follows that 


( 11 ) 


Q{k,m,a,b,p) = lini g{x) 

x-^-oo 


dx 

xePo= 


lim g(x) 

x—>0 


x^ dx r 

xeP- Jo [J 


xk 1 
gpx 



—Qm(a, 6vx)da: 


( 12 ) 


where 


g{x) = Qm{a,by/x) 


(13) 


The integrals in dTB ean be expressed in terms of the ineomplete gamma funetion yielding 
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g{k,m,a,b,p) = lim - lim 

x^oo p'^ a:^0 p^ 

Recalling the identities 


1 r°° 

T / mpx) 


—Qm{a,by/x)dx. 

ax 

(14) 


r(a, oo) = 0 , 


(15) 


r(a,0) = r(a) (16) 

and 

Qm{a,0) = l (17) 

as well as setting u = y/x in (fT4l) along with utilizing ^ eq. (9)] and [l40l eq. (8.352.4)], it 
follows that 


g{k,m,a,b,p) = 


T{k) 


X 


pK 


k-l 


1 - Y, ^ I x^-^^e-(^+-^yim-iiax) dx 


1=0 


(18) 


The above integral can be expressed in closed-form using HOl eq. (6.621.1)]. To this effect, eq. 
([TS]) is deduced, which completes the proof. ■ 


B. Closed-form Solutions to lF{k, m, a, b, p) 

A closed-form expression for (|2|) was reported in [|2T]| for the case that k is integer and m is 
arbitrary real, namely 


J^(/c, m, a, 6,p) 


r(i)r (m, f) 

p^r(m) 

^2^2mr(fc)e-f ^ / 2p \' 
m\p^2'^{a? -f 2p) \a^ 2p) 


iFi 


-f 1; m -f 1; 


a^b^ \ 
2a? -\- Ap) ' 


(19) 


Likewise, in what follows we derive a closed-form expression for the useful case that m is 
integer and k is arbitrary real. 
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Theorem 3. For k,p E M'*', m G N and a,b eM, the following closed-form expression is valid 


r(i-) 2‘r(i)e-^<i.Ji,i,i;f.5#Sj) 

J^{k,m,a,h,p) = — - 


pft, 


(a2 + 2pf 


( 20 ) 


r(fc)2^e-^ / ^2^2 


[a? + 2p)^ n!2 


20^ + 4p 


where $ 2 (') is the Humbert function of the second kind 


Proof: The series in |[TS eq. (9)] can be also expressed as 


J-{k, m, a, b,p) 


n=o n!b~'^’^2’^~^(a^ + 2p)^e^ 

r(fc)2^e-^ b^^ f a%^ \ 

(a^ + 2pY n\2'^ ^ ^ V + 4p J 

"- — -V-^ 

'H{k,a,b,p) 


Applying [l40l eq. (9.14.1)] in 'H{k,a,b,p) and given that 


( 21 ) 


one obtains 


(n + 1)/ = 


r(n + / + 1) 
r(n + 1) 


(l)n+Z 

(l)n 


( 22 ) 


H{k, a, b,p) 


e-^2"r(fc) il)nik)i ^ ( 2 ^^) 

(a2 + 2p)fc (l)n+z n\ l\ 


(23) 


Notably, the above expression can be expressed in terms of the Humbert hypergeometric function 
of the second kind yielding 


FL{k, a, b,p) 


e-"^2^T{k) (. j a^b‘^ \ 

(a2 + 2p)^ V^’ Y’ 2a2 + 4p ) ‘ 


Inserting (l24l) into (l2TI) yields (l20l) . concluding the proof. 


( 24 ) 


C. Specific Cases of lF{k,m,a,b,p) and Q{k,m, a,b,p) 
Simple expressions are derived for /c = 1, a = 0 and 6 = 0. 
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1) The case that k = 1: In this special case it follows that 


a,6,p) = / Qm{a,h^/x)e 

Jo 


(25) 


and 


J^(l, m, a, 6, p) = / Qm{(i\/xib)e ^^dx. 


(26) 


Equation (1251) is given by ifT^ eq. (16)]. In the same context, a generic closed-form expression 
for (| 2 ^ is derived below. 

Lemma 1. For a, m G M and b,p ^ M'*', the following closed-form expression is vali^ 


J^(l, m, a, b,p) = 


r(m.¥ 


+ 




pb^ 

,2-zjt 


pT{m) pa^"*r(m)(a^ -f 2p)^" 
where 7 ( 0 , x) denotes the lower incomplete gamma function. 

Proof: By setting A; = 1 in (fT9l) . it immediately follows that 


J^(l, m, a, 6,p) = 




pr(mj m!p2™'(a^ + 2p)e"^ 

Notably, using the following hypergeometric function identity 


(27) 


(28) 


iFi(l; n; x) = "(n — l)e'^ 7 (n — 1, x) 


(29) 


and recalling that 


7 ( 0 , x) = r(a) — r(a,x) 


(30) 


and 


Tfl ^ 

—^ = (m — 1)! = r(m) (31) 

m 


^The integral in 1201 eq. (1)] was recently evaluated in closed-form. However, this solution does not account for m, a, b, p) 
since the two integrals would be equal only when /o(0) = 1. Yet, this can be achieved for p 2 ~ t and c = 0, which eliminates 
the power term and as a consequence, the integral reduces to ( I26t . which is simply expressed in closed-form in t27l (. 










equation ( 1281 ) can be equivalently expressed as 


J^(l, m, a, b,p) = 


T j (a^ + 2p)"*“^e ^+2p (a^ + 2p)'^~^e “'^+2? / 

pr(m) pa‘^^~‘^ 2a^ + ^ 


20^ + 4p 

(32) 


To this effect and performing long but basic algebraic manipulations (|3^ reduces to (1271) . which 
completes the proof. ■ 

Remark 1. A similar expression can be obtained through (l20l) . 

2) The case that a = O.' In this special case it follows that 


and 


Q{k,m,0,b,p) = / ^Qm{,^,b^/x)e 

Jo 


J^{k,m,0,b,p) = Qm{d,b) / x^ ^^dx. 

Jo 


(33) 


(34) 


Lemma 2. For k,p E M"*", m eN and 6 G M, the following closed-form representations hold 


^(fc,m, 0 , 6 ,p) = 


m —1 


b‘^T{k + 1) 
(JF + 2pY + 2pY 


(35) 


and 


r(A;) 


J^(fc, m, 0 , 6 ,p) = 

Proof: 

Applying eq. (2)] in (1331) one obtains 


-e 2 


pK 


/! 2 '' 


1=0 


(36) 


^(/c,m, 0 , 6 ,p) = ^ 


1=0 


l\2^ 


X 


k+i-i -x{p+^) 


dx. 


(37) 


Evidently, the above integral can be expressed in closed-form in terms of the Euler gamma 
function. This is also the case for J^(/c, m, 0, &,p) since Qm{d,b) is not a part of the integrand. 
As a result, by performing a necessary change of variables in [l40l eq. (8.310.1)] and substituting 
(1371) and (l33l) one obtains (1351) and (l34l) . respectively, which completes the proof. ■ 











3) The case that b = 0; In this special case it follows that 


'OO 


Q { k , m , a , 0 , p ) = Qm { a , 0 ) 


^k-ie-px^^ 


(38) 


and 


'OO 


J ^{ k , m , a , 0 , p ) = 


X 


.k-l 




(39) 


Lemma 3. For p G M’*' and a , m,k E M, the following simple closed-form expression is valid 



(40) 


Proof: The proof follows with the aid of the identity Qm(a, 0) = 1, in [|3 eq. (1)] as well 


as EOl eq. (8.310.1)]. 


III. Applications in Wireless Communications 


As already mentioned, the derived expressions for ([I]) and ^ can be used in applications 
relating to natural sciences and engineering, including wireless communications and signal 
processing. Based on this, they are employed in the derivation of simple expressions for the 
detection of unknown signals in cognitive radio and radar systems as well as for the channel 
capacity using switched diversity. To this end, a closed-form expression is firstly derived for 
the average probability of detection over Nakagami—m fading channels. Likewise, a closed- 
form expression is derived for the channel capacity with channel inversion and fixed rate in 
switch-and-stay combining (SSC) under correlated Nakagami—m fading conditions. 

A. Energy detection over Nakagami— m fading channels with arbitrary values of m 

The detection of unknown signals is modeled as a binary hypothesis-testing problem, where 
FIq and Hi denote the cases that a signal is absent or present, respectively. The corresponding 
test statistic is typically represented by the central chi-square and the non-central chi-square 
distributions, respectively, and is compared with an energy threshold, A [HI. 
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Fig. 1. ROC curve for u = 5, 7 = 15dB and different values of m. 


Corollary 1. For 7 , A G M’*', and either m > 0.5 and u G N, or m ^ N and u G M"*", the 
average probability of detection over Nakagami—m fading channels can be expressed as 

Pd = . P im,u,V2,'/X,— ]. (41) 

7 ^r(m) \ 7 / 

Proof: The probability of false alarm and probability of deteetion in additive white Gaussian 
noise are given by Pf = r(M, A/2) and Pd = Qu{s/^^ V^), respeetively, where u and 7 denote 
the time-bandwidth produet and the instantaneous signal-to-noise ratio (SNR), respeetively [[20ll . 
ETI . If2^ . It is reealled that in energy deteetion over fading ehannels, the Pd is averaged over 
the fading statistics. To this effect, for the case Nakagami—m fading channels in I® eq. (2. 21)] 
is represented as follows 

_ ™rn, poG _ ^ 

Pd = (42) 

7™r(m) Jo 

Evidently, the above integral can be expressed in terms of (fT9l) or (l20l) . This yields (HTI) . which 
completes the proof. ■ 
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Notably, the offered expression can account for arbitrary values of m, contrary to existing 
analyses that assume integer values of m for simplicity. Fig. 1 illustrates the corresponding 
probability of missed detection versus probability of false alarm ROC curve for different values 
of m. One can notice the sensitivity of m, particularly for small values, and thus, the usefulness 
of the offered expression also in practical scenarios. 

B. Capacity with channel inversion and fixed rate over correlated Nakagami—m fading using 
switch-and-stay combining 

Channel capacity under different transmission policies is particularly useful in achieving certain 
quality of service requirements. In this context, channel inversion with fixed rate (CIFR) has 
been rather useful as it ensures a constant SNR at the receiver through adaptation of the transmit 
power. This method relies on fixed-rate modulation and fixed code design, which renders its 
implementation relatively simple HI . 



Fig. 2. Average SE vs 7 for CIFR under correlated Nakagami—m fading with SSC for 7 t = OdB, p = 0.5 and different 
values of m. 
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Corollary 2. For 7 , -B G M"*", m > cind 0 < p < 1, the capacity with channel inversion and 
fixed rate over correlated Nakagami—m fading channels with SSC can be expressed as 


where 


CciFR — B log 2 



n{m,y,yT,p) 


(43) 


2myT 

(1 -p)7’ 

(44) 

with p and denoting the correlation coefficient and the predetermined SNR switching thresh¬ 
old, respectively. 

Proof: The CIFR over fading ehannels is defined as iHTI 

CciFR = B log 2 1^1 + j . (45) 

In the ease of switehed diversity and eorrelated Nakagami—m fading, the PDF of the SSC output 
is given by I® eq. (9.334)]. By also setting 

n= 

Jo 7 

it follows that 


n{m,y,yT,p) = 


m 


yV{m) 


m 


1 , 


myT 


m 


7 / y^T{m) 


m — 1, m, 


2mp 

(1 -p)7’ 


11= A{y)dy+ / p-y{y)dy 
•Jo •J '7'j' 


(47) 


where ^ 4 ( 7 ) is given in ||4l eq. (9.335)]. To this effeet and using flU eq. (2. 21)] one obtains 


^("i,7,7T,p) 


m 


^r(m) 


7 ™ T" d7 


'It 


- r 1 d7. 


(48) 


/o T"r(m) \^y (1 - p)7’ Y (1 - p)7^ 

The first two integrals in (148]) ean be expressed in terms of the gamma funetions, whereas the 
third integral has the algebraie form of dH). Therefore, by performing the neeessary ehange of 
variables yields (l43l) . whieh eompletes the proof. ■ 


^l| 3 
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The behavior of the eorresponding average speetral effieiency versus average SNR is illustrated 
in Fig. 2 for different values of m with fixed values of 7 t and p. The signifieant effeet of the 
severity of fading on Ccifr is elearly observed. 

IV. Conclusion 

Novel closed-form expressions were derived for two Marcum Q—function integrals which are 
both simple and generic. Simple analytic expressions for involved special cases were also derived 
in closed-form. These expressions are tractable and are expected to be useful in analyses relating 
to natural sciences and engineering, including wireless communications and signal processing. 
To this end, they were employed in the analysis of energy detection in RADAR and cognitive 
radio systems as well as in the channel capacity with channel inversion and fixed rate over 
correlated multipath fading channels. 
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